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We investigate the fine structure of graphene on iridium, which is a model for graphene weakly 
interacting with a transition metal substrate. Even the highest quality epitaxial graphene displays 
tiny imperfections, i.e. small biaxial strains, ca. 0.3%, rotations, ca. 0.5°, and shears over distances 
of ca. 100 nm, and is found incommensurate, as revealed by X-ray diffraction and scanning tunneling 
microscopy. These structural variations are mostly induced by the increase of the lattice parameter 
mismatch when cooling down the sample from the graphene preparation temperature to the mea- 
surement temperature. Although graphene weakly interacts with iridium, its thermal expansion is 
found positive, contrary to free-standing graphene. The structure of graphene and its variations are 
very sensitive to the preparation conditions. All these effects are consistent with initial growth and 
subsequent pining of graphene at steps. 



I. INTRODUCTION 

Graphene preparation at the surface of low-carbon 
solubility metals like Irji Cuj 2 - or Pt£ is a surface- 
confined process which stops once the surface is pas- 
sivated by a full graphene layer. Therefore it is 
a straightforward route towards the production of 
large-area, highly quality graphene. Not only is the 
carbon solubility low in the aforementioned metals, 
but also their interaction with graphene. In particu- 
lar, in the model graphene on Ir(lll) system, where 
graphene was shown to be almost free-standing, 4 
some of the longstanding issues in graphene re- 
search could recently be addressed, noteworthy the 
engineering of graphene's Dirac cone via bandgap 
opening-^ or Fermi velocity renormalization^ elec- 
tron confinement in graphene quantum dots^r— or 
the use of graphene as a weakly-perturbating and 
protective layer for fragile surface phenomena such 
as the Rashba-split Ir surface stated The weak in- 
teraction of graphene with metals such as Ir, Cu or 
Pt also has important consequences regarding the 
structure of graphene. Indeed, graphene domains 
with different stacking with respect to the substrate 
tend to form on these metals, presumably because 
the formation of each of the domain involves sim- 
ilar energetic costs. This was first observed some 
decades ago^i 2 - further addressed in details on Ir ^ 13 ! 14 
and latter on, on Cu^ and Pt At the bound- 
ary between these domains, dislocations (heptagon- 
pentagon pairs) are found. 13 

Such structural defects hinder the improvement 
of graphenes performances, e.g. for electronic 
transport ^ 19 ! 20 towards state-of-the-art ones ob- 



tained for exfoliated, suspended graphene. Cer- 
tain defects are on the contrary desirable: some 
point defects are expected to switch on magnetism in 
graphene^i inhomogeneous strain fields were shown 
to induce electronic gap opening, 22 and defects are 
thought a pathway for the intercalation of various 
species, which recently appeared a powerful route for 
preparing graphene hybrid systems^ 3 - — Producing 
defect-free single-crystalline graphene or graphene 
with well-known defects having specific properties on 
metals such as Cu or Ir is thus of both fundamental 
and applied interest. Recently, strategies have been 
developed in order to avoid the formation of grain 
boundaries, in graphene grown on Ir(lll) i 27 i 28 Scan- 
ning tunneling microscopy (STM), low energy elec- 
tron diffraction, and low-energy electron microscopy 
allowed concluding to a single crystalline orientation 
in graphene at the scale of a centimeter; however, 
the question of graphene's crystalline perfection can 
only be answered with higher resolution probes, such 
as X-ray scattering, which is also an unvaluable tool 
to gain knowledge about the nature of defects. 

We report here the investigation of the structural 
properties in graphene/Ir(lll), by in situ high reso- 
lution X-ray diffraction, which allowed us to unveil, 
together with STM, small deviations to a perfect 
graphene lattice and provide insight into the nature 
of the epitaxy between graphene and the metal. 

The article is organized as follows. After present- 
ing the experimental methods, we turn to the experi- 
mental results, and successively address the average 
lattice parameter in graphene, the distributions of 
the structural parameters of graphene, and the de- 
pendence of these parameters on temperature. The 



last section of the article is devoted to the discussion 
of these results. 



II. METHODS 

X-ray diffraction with synchrotron light was con- 
ducted inside the ultra-high vacuum growth cham- 
ber, under grazing incidence to achieve maximum 
sensitivity to the graphene overlayer; STM was per- 
formed in a separate system. 

The iridium single crystal polished in a (111) 
plane, bought from the "Surface Preparation Lab- 
oratory", was prepared by repeated cycles of 0.8- 
1.5 keV Ar + sputtering at room temperature (RT) 
and high temperature flash. In order to decrease the 
crystal mosaic spread and increase the terrace width 
(of the order of 100 nm after surface preparation), 
flash temperatures as high as 1770 K were employed. 
Repeated sample annealing at 1070 K under O2 par- 
tial pressures of 10 -8 mbar, during 2 h, allowed re- 
moving carbon incorporated into the sample during 
high temperature flashes. The sample temperature 
T s was measured with a pyrometer (60 K uncertain- 
ties). 

Graphene growth was performed by first adsorb- 
ing ethylene at RT during 5 min with a 10 _T mbar 
pressure of ethylene in the chamber, then flash- 
ing the sample temperature to 1470 K, and finally 
decreasing the sample temperature to T g = 1070- 
1270 K and letting in a 10 ~ 8 mbar pressure of ethy- 
lene for 10 minj2£ This is known to yield a single 
layer graphene sheet with almost 100% coverage, 
as shown on the atomic force microscopy images in 
Fig. [TJ free of defects such as grain boundaries, with 
a single crystallographic orientation, at least to the 
accuracy of STM and LEED. Compared to growth 
above 1270 K without the first step (RT absorption 
and flash anneal at 1470 K), it also reduces the lat- 
tice mismatch between graphene and its substrate 
at RT after growth, and thus limits the density of 
wrinkles which form upon cooling. Platinum was de- 
posited at RT on graphene/Ir(lll) using electron- 
beam evaporation of a high purity Pt rod under 
ultra-high vacuum. 

Grazing incidence X-ray scattering (GIXS) was 
performed at the INS instrument installed at the 
BM32 beam-line of the European Synchrotron Ra- 
diation Facility (Grenoble, France), in a ultra-high 
vacuum system coupled to a Z-axis diffractometer 
and with a base pressure below 10 -10 mbar. Several 
experiments were performed with a monochromatic 
photon beam of either 11 and 21 keV incident under 
angles of 0.38 and 0.19° respectively. The corre- 




FIG. 1. (Color online) (a,b) Atomic force microscopy im- 
ages measured in atmospheric conditions on an Ir(lll) 
single-crystal covered with graphene grown under ultra- 
high vacuum. Note the interconnected network of wrin- 
kles, which reveals the presence of graphene over the 
whole surface. 



sponding X-ray attenuation length for Ir is 57 A at 
11 keV, and 68 A at 21 keV. The incident beam was 
doubly focused to a size of 0.4x0.3 mm 2 (full width 
at half-maximum in horizontal and vertical direc- 
tions, respectively) at the sample location. Detector 
slits, located 640 mm away from the sample, were set 
at 1 mm (or 0.5 mm) parallel to the sample surface 
(which was vertical) and 8 mm perpendicular to it, 
resulting in angular acceptances of 0.09° (or 0.045°) 
in the vertical direction and 0.7° in the horizontal 
direction. GIXS data were measured with the help 
of a standard Nal scintillation point detector. 

When used, the reciprocal lattice units (r.l.u.) ft, 
fc, and I are calculated using a hexagonal surface 
unit cell with lattice vectors ai 5S = l/2(ai jV + a2 jV )> 
a 2 , s = l/2(ai )V -a 2 , v ), a 3)S = ai )V +a 2 , v +a 3)V , where 
ai, v , &2,v a3 ' v are ^ ne v °l ume urn t cell lattice vectors 
of Ir, all having 3.839 A length, so that surface lat- 
tice vectors have 2.714 A length in-plane and 6.649 A 
out-of-plane. All scans of the scattered intensity 
which are plotted and discussed hereafter have been 
acquired at an out-of-plane reciprocal space lattice 
unit I = 0.09. 

The detector slits (opened 1 mm, located at 
640 nm from the sample), intercept a 2tt/X x 1/640 
region in reciprocal space along the slits direction, 
A being the wavelength of the X-ray beam. They 
are inclined by an angle which equals the scattering 
angle with respect to the radial scan direction (uj). 
Due to the slits, the peaks are broadened in radial 
and azimuthal directions. This broadening varies as 
a function of the peak order and is systematically 
corrected for in the following. 



2 



III. AVERAGE LATTICE PARAMETER IN 
GRAPHENE 



IV. DISTRIBUTION OF LATTICE 
PARAMETER AND CRYSTALLOGRAPHIC 
ORIENTATION IN GRAPHENE 



The in-plane structure of graphene was investi- 
gated by exploring an in-plane cut in reciprocal 
space (Fig. [2j), i.e. measuring the scattering inten- 
sity as a function of the radial Q r component of the 
scattering vector and of the azimuthal angle uj. In 
order to test the structural quality of graphene, spe- 
cial attention was paid to the moire resulting from 
the difference in lattice parameter between graphene 
and Ir(lll), which is known to amplify tiny deforma- 
tions and rotations^ and is associated to diffraction 
peaks as we will see. 

At locations where Ir crystal truncation rods 
(CTRs) or graphene diffraction rods intersect the in- 
plane cut of the reciprocal space, intensity maxima 
are found. For the Ir CTR passing through the cen- 
ter of the second (Tj r ) and third (r 2 r ) Brillouin zones 
(BZs), the result is shown in Fig. [3j The Ir(lll) sin- 
gle crystal yields sharp contributions reflecting its 
high quality. Besides the Ir peaks, a graphene con- 
tribution is also observed (at and T^), whose 
Q r position yields the RT in-plane projection of the 
lattice parameter of graphene, according to ay = 
47r/(v / 3Q r ), 2.4543±0.0005A. 

A peak, labeled in Fig. [3^ and referred to as 
a moire peak in the following, is found halfway be- 
tween the Ir (r 2 r ) and graphene (r^) peaks, cor- 
responding to a a m = 25.6±0.2A period. The ori- 
gin of this peak is discussed in section IVII Higher 
order moire peaks, which are revealed by electron 
diffraction 29 thanks to the strong interaction of elec- 
trons with matter giving extreme surface sensitiv- 
ity and multiple diffraction, are not measurable here 
around r 2 r ; however, as sketched in Fig. [2j other first 
order moire peaks were measured around Fj r c and 

r 4 

1 Ir,C 

Platinum deposition on graphene/Ir(lll), which 
yields nanoclusters self-organized on the moire^ in- 
creases the intensity of the moire peaks and makes 
new ones appear (Fig. ^jp). We note that in the 
second BZ, the Tq position is also the position of 
the first order moire peak around Tj r . This is pre- 
sumably the origin of the increased intensity at this 
location in the presence of Pt clusters. The presence 
of Pt does not alter the position of the graphene 
peaks (Fig. [3}d). 



The line-shape of diffraction peaks is sensitive 
to crystal imperfections i^I finite size effects, mosaic 
spread (in the plane of the sample surface consid- 
ering the diffraction geometry which we employ), 
and distributions of lattice parameters all yield peak 
broadening, the origin of which can be discrimi- 
nated by studying the broadening dependence on 
the diffraction peak order. Along the radial direc- 
tion, the graphene peak width is found to increase 
linearly with peak order from Tq to Tq (Fig. |4^i), in 
contrast with the Tj^ 4 Ir peaks, whose radial width 
is constant (0.01 A -1 ) with order. This implies that 
the lattice parameter in graphene is spatially vary- 
ing, with a full width at half maximum (FWHM), 
given by the increase of the peak width with of peak 
order (Fig. |4]d), 0.0093±0.003 A. The FWHM also 
has a constant component, obtained by linear ex- 
trapolation to zero order (center of the reciprocal 
space), 0.011 A -1 , which corresponds to a finite size 
effect pointing to structurally coherent graphene do- 
mains of 60±30 nm extension. 

While their radial width increases as a function 
of peak order, the graphene peaks have constant 
azimuthal width (0.53°, substantially larger than 
the constant 0.08° azimuthal width of Ir peaks), as 
sketched in Fig. [2] and shown in Fig. [4k, which in- 
dicates that small- angle disoriented and/or sheared 
domains are also present. Similar rotations were 
reported recently in graphene on Ru(0001), using 
low energy electron microscopy and micro-LEED^ 2 - 
Such rotations should translate into a ca. 10 times 
larger rotation of the moire peaki^ with respect to 
the nearest Ir peak, i.e. into an angular width of the 
moire peak with respect to the origin equal to that 
of the graphene peaks divided by its order (Fig. [2^). 
The measured 0.26° FWHM of the moire peak 
confirms this interpretation. Note that the coex- 
istence of in-plane rotated phases was already re- 
ported in the case of large-angle orientation variants 
in graphene/Ir, 14 and small- angle scatter between 
graphene islands on Ir^ 

In the present study, an attempt to further de- 
tect lattice variations or distortions by STM could 
only reveal shears. Shears were detected in STM 
by scrutinizing the alignment and misalignment on 
carbon atomic rows and of moire high symmetry di- 
rections. Figure [5] highlights a triangular region of 
graphene/Ir (111), into which the carbon zigzag rows 
align to the high symmetry directions of the moire, 
only on one of the sides of the regions. The mis- 
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FIG. 2. (Color online) Sketch of the first to fifth graphene BZs and their centers for Ir (T^ 4 ) (red) and graphene 
(r^ 4 ) (blue) where the scattered intensity was measured. Blue ellipses represent in-plane cuts of graphene diffraction 
rods, characterized by their position, yielding the average projected lattice parameter, as well as by their Q r radial 
and uj azimuthal widths. For the sake of clarity, the moire diffraction peaks have been omitted from this sketch, 
but are shown in the close-up view, which marks the measurable moire peaks and gives a geometric construction 
explaining the observed width of these peaks. 
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19 Ir* 21 C '"9 Ir* 10 C ft(r.ku.) 

FIG. 3. (Color online) (a) In-plane radial scans at RT close to the third BZ center for Ir (r^) and graphene (Tq) 
revealing the Ir CTR (red), the graphene rod (blue) and a moire peak (T^, gray), as a function of radial momentum 
transfer (Q r ) and reciprocal space lattice units (h). The position of the commensurate structures with 9 Ir cells 
matching 10 C ones and 19 Ir cells matching 21 C ones are marked, (b) In-plane radial scans measured for another 
sample (thus having different structure, see section E}, with and without Pt clusters on top (square and disk symbols, 
respectively), close to the second order BZ center for Ir (Tj r ) and graphene (r^), revealing that moire peaks (T^) 
appear and that the graphene peak is reinforced due to the presence of Pt clusters. Peaks in (b) include two 
components originating from two Ir crystallites. Tq is thus fit with two gaussians, whose centers do not change upon 
Pt cluster deposition. Insets: sketch of the moire reciprocal lattice. 



alignment on the two other sides is of the order of a 
degree. It is typically one order of magnitude smaller 
between the graphene and Ir(lll) atomic rows. This 
small shear is in the range estimated by measur- 
ing the azimuthal width of the diffraction peaks of 



graphene. As we will explain in section fVTl detecting 
the distribution of lattice parameter in graphene by 
STM is an extremely demanding task, beyond the 
scope of this article. 




FIG. 4. (Color online) (a) Radial scans for graphene at different peak orders (Tq~^ 4 ), vertically shifted for clarity, (b) 
Graphene radial scans FWHMs (AQ r ) as a function of peak order (T^ 4 ), at RT and 1070 K, and linear fits (lines), 
(c) Azimuthal angle scans at different peak orders (rj^ 4 ), vertically shifted for clarity. For T^, the azimuthal angle 
scan (blue disks) has been smoothened (solid lighter blue line) for easing the assessment of the scan FWHM. (d) 
Graphene radial (AQ r ) and azimuthal (Aoj) FWHMs as a function of the sample temperature for the second order 
peaks. 



V. DEPENDENCE OF GRAPHENE 
STRUCTURE ON TEMPERATURE 

A. Dependence on growth temperature 

The average epitaxial matching of graphene on 
Ir(lll) was found to depend sensitively on T g . Fig- 
ure [6^ shows a graphene peak (light blue) whose 
position corresponds to a||=2.4470±0.0005 A, 0.3% 
smaller than the previous value, obtained with an- 
other sample (see Fig. [3^)- The accuracy of our tem- 
perature probe (±60 K) does not allow detecting dif- 
ference in T g between the two samples. Actually, the 
RT graphene lattice parameter a\\ was found to vary 
by almost 0.01 A for different preparations, with no 
clear dependence on T g (Fig. [6}d). 

Discrepancies were found also in the value of the 
radial width of the graphene peaks, between two 
samples tentatively prepared at the same tempera- 
ture, within the accuracy of the temperature probe. 
These variations are accounted for in Fig. 0]i in the 
form of error bars. The general trends, namely the 
FWHM slope and extrapolation to the center of re- 
ciprocal space, do not vary significantly. This im- 
plies that the two samples have similar domain size 
and width of the distribution of graphene lattice 
parameters. We found no clear dependence of the 
graphene radial width neither on T g nor on ay, in 
line with the behavior of ay versus T g . 



B. Dependence on sample temperature 

Varying T s away from T g also influences the epi- 
taxial matching: while between RT and 1070±60 K 
the Ir lattice parameter is expanded by 0.48%, 
for graphene a\\ is increased by only 0.33%, from 
2.4470±0.0005 to 2.455±0.001 A (Figs. Ek,b). The 
superperiodicity is slightly decreased, to 24.5±0.2 A. 
The radial scans widths (Fig. HJd) are again found to 
increase with peak order, but to a much lesser extent 
than at RT. The corresponding distribution of lattice 
parameter, 0.0023 A, is almost one order of magni- 
tude smaller than at RT, but the extrapolation to 
zero of the radial scans width yields a value for the 
size of the structurally coherent graphene domains 
similar to the RT value, 60±30 nm. The angular 
widths are also constant (0.43°) as a function of or- 
der, but smaller than at RT. Figure 2]i shows that 
both the radial and angular widths of a graphene 
peak increase with decreasing the sample tempera- 
ture. 



VI. DISCUSSION 

The RT lattice parameter of isolated graphene 
-a conceptual object- was calculated to be <2o,rt 
= 2.4565 A^ In the above we have extracted 
the in-plane projection of the lattice parameter in 
graphene, ay. Deducing from this parameter the 
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FIG. 5. STM topograph of graphene/Ir(lll) measured 
with a 10.9 nA and a 0.2 V tunneling current and bias 
respectively. Solid lines follow the moire lattice, white 
dots mark centers of carbon rings. The two triangles 
share only one edge, evidencing an in-plane shear. 



actual lattice parameter, <2rt, requires taking the 
effect of the graphene/Ir moire into account. The 
weak, varying interaction between graphene and 
Ir as a function of the location in the moire (see 
Fig. [7)) is indeed associated with a nanorippling 
of graphene, characterized with a 5=0.5 A typical 
height amplitude A Assuming a sinusoidal variation 
of the height, <2rt is typically 0.3% larger than ay. 
This yields a c =2.4617A and 2.4548 A for the sam- 
ples studied in Figs. [3^ and [6^ respectively, i.e. 0.2% 
larger and 0.1% lower than ao,RT respectively. The 
self-organization of Pt clusters on top of graphene 
does not modify the in-plane component of the lat- 
tice parameter of graphene, however it could affect 
the nanorippling, in turn changing its lattice param- 
eter. The small deviations of ac from <2o,rt are con- 
sistent with the weak interaction between graphene 
and Ir(lll), characterized by a limited charge trans- 
fer between graphene and Ir, typically 0.01 electron 
or hole per atom. 4 

By contrast, in graphene/Ru(0001), a much 
larger ay =2.4895 A value was founds In this sys- 
tem, nanorippling is more pronounced than in 
graphene/Ir (111), corresponding to an expected 
height variation of 1.5 A—, which yields ac roughly 



2% larger than &o,rt? consistent with the strong elec- 
tron dopin g 36 ! 37 of graphene on Ru. 

The observed weak moire peaks are associated 
to the weak but non vanishing graphene-Ir(lll) in- 
teraction. We assume a small in-plane sinusoidal 
strain in Ir (sketched in Fig. [7] with a color gradient 
for Ir atoms) stemming from this interaction, hav- 
ing locally a weak tendency to covalent bonding, 4 
which according to density functional theory calcu- 
lations moves Ir atoms away from their position in 
bare Ir(lll) by a e value of a few 0.001 A with a 
a m period 3 ^ as a function of their position in the 
moire lattice. In a simple one-dimensional model, 
from the position of an Ir atom labeled i in a chain, 
ai r x (i + £COs(27riai r /a m )), the scattered amplitude 
calculated by summing up the contribution of all 
Ir atoms accounts for the observed so-called moire 
peaks, whatever the relation, commensurate or not 
between the Ir, graphene and moire lattice param- 
eters. The periodic modulation of ac in graphene 
also contributes to the moire peak intensities, but to 
a much smaller extent, since scattering by C atoms 
is ca. 40 times lower than by Ir ones. 

The measured ay values are substantially larger 
than the value (2.4435 A) expected for a com- 
mensurate phase with a moire consisting of 10 
graphene cells matching 9 Ir ones. Second 
order commensurability, with 25 graphene unit 
cells matching 23 metal ones, was concluded for 
graphene/Ru(0001)^ The closest high order (sec- 
ond and third) graphene/Ir (111) commensurate 
moire would consist of 21 graphene unit cells match- 
ing 19 metal ones, which is yet slightly but signifi- 
cantly off the a m and ay values. We conclude that 
on the average, at RT, the graphene and Ir lattices 
are incommensurate. 

However, a description in terms of a single in- 
commensurate graphene/Ir (111) phase would be too 
naive as a distribution of lattice parameters was un- 
veiled. This distribution is plotted in Fig. [TJ It is 
tempting to invoke the periodic modulation of ay due 
to the moire nanorippling (see Fig. [7]) as the interpre- 
tation of this distribution. Nevertheless, these varia- 
tions are periodic (or quasi-periodic in case of incom- 
mensurate structure) and thus contribute to all ob- 
served peaks without yielding any broadening of the 
graphene peaks. Small random or uncorrelated dis- 
placements similar to thermal motion, either static 
or dynamic, are also excluded as they only yield a 
Debye- Waller like decrease of the peak height with 
order, without peak widening^ 1 - Only two scenarios 
are actually relevant for accounting for the distribu- 
tion of lattice parameter: the existence of domains 
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FIG. 6. (a) Radial scans close to the third BZ center at 1070 K (light red and light blue), and the same sample at 
RT (red and blue). Note that the RT measurement corresponds to a sample decorated with Pt clusters reducing 
counting times (see Fig. 02)). (c) Spread of in-plane projections, ay, of ac, measured at RT, as a function of the 
growth temperature T g , for two series of measurements (yellow and red data points). The orange point corresponds 
to a sample grown at 1170 K and studied at 1070 K. 



with different average lattice parameters (first sce- 
nario) or progressive changes of lattice parameter 
within large domains (second scenario). These sce- 
narios are discussed more into details at the end of 
this section in the light of the effects of temperature, 
which we will now address. 

The fact that the RT lattice parameter ac varies 
by several tenth of a percent as a function of the 
growth temperature T g , without simple relationship 
with it, suggests that due to the mismatch in thermal 
expansion coefficient (TEC) between graphene and 
Ir^ the epitaxy of graphene on Ir(lll) is initially 
set at T g , and that the strain induced during cooling 
down to RT by the mismatch in TEC might modifies 
the epitaxy in a non trivial way. 

Taking nanorippling into account as before, we 
find that at 1070±60 K ac is 0.3% larger than 
the expected lattice parameter (ao,io70 = 2.4556 A) 
for isolated graphene at this temperature^ At RT 
for this sample, the deviation to the expected <2rt 
value was less than 0.1%. The average epitaxial 
relationship thus indeed changed between 1070 K 
and RT, which suggests some slippage of graphene 
on Ir(lll) against temperature-induced stress, pre- 
sumably due to the high stiffness of graphene. The 
TEC of graphene is found positive, at variance 
with the vanishing 34 or negative 4 ^ lattice param- 
eter variation predicted for free-standing graphene 
in this temperature range. The TEC is much 
larger than estimated from first principle calcu- 
lations for graphene/Ir(lll)j 4 i The partial inheri- 



tance, in graphene, of the Ir TEC, could result from 
the weak tendency to periodic chemical bonding in 
graphene/Ir (Ref. 0), or from more subtle effects 
near steps edges (see below). A more extreme situ- 
ation was encountered in a related system, fo-BN on 
Rh(lll) (Ref. EH) in which the two TECs were found 
equal presumably due to the very strong interaction 
between the two materials. 

Graphene not only exhibits domains with differ- 
ent structures (see discussion below) at RT, but also 
at 1070 K (Fig. 0b). The fact that at this tempera- 
ture wrinkles are absent 39 suggests that the domains 
develop during growth. The spread in the structure 
between domains is however amplified as the tem- 
perature decreases from T g (Fig. |4ji). This suggests 
that graphene minimizes its elastic energy by am- 
plifying the structural variations as the compressive 
epitaxial stress, due to a smaller TEC than that of 
Ir, increases. 

This brings us to the discussion about the two 
above-mentioned possible scenarios. 

The first scenario is sketched in Fig. [8^. It shows 
domains with a size of ca. 60 nm, most being in- 
commensurate with Ir(lll), possibly coexisting with 
commensurate phases (21 C hexagons coinciding 
with 19 Ir distances for sample of Figs. [3bt and 11 
C hexagons matching 10 Ir for sample in Fig. [6^) 
present in small proportion. The coexistence of var- 
ious graphene phases in this scenario implies the ex- 
istence of a large number of shallow minima in the 
energy of graphene upon nanorippling, bond strech- 
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FIG. 7. (Color online) Construction of the distribution of 
projections of lattice parameters in graphene, {ay } , from 
the distribution of lattice parameters, {ac}, taking the 
nanorippling of graphene along the moire into account. 
The purple histogram shows a distribution of graphene 
lattice parameters randomly chosen in a standard nor- 
mal (Gaussian) distribution (as a guide for the eye, the 
corresponding gaussian function is plotted with a purple 
line), with FWHM corresponding to the increase of the 
width of the graphene peak as a function of BZ order in 
the radial direction (Fig. \4]p). The grey histogram is the 
distribution of in-plane projections yielded by nanorip- 
pling for one specific value of the lattice parameter. The 
cyan histogram includes all such distributions each cor- 
responding to all lattice parameter values included in the 
distribution shown in the purple histogram. The distri- 
butions were derived by taking into account 10 moire 
unit cells in a one-dimensional model. The inset shows a 
side-view of this model along a (112) direction in Ir illus- 
trating the nanorippling- induced {a^i} distribution for 
a specific ac , and the periodic strain in the last Ir plane 
(red shades, s amplitude). The characteristic lengths of 
the model, the moire lattice parameter a m , the graphene 
lattice parameter a, its in-plane projection a\\, and the 
nanorippling amplitude s are shown. The V3 factor cor- 
responds to the (112) direction. 



ing/compression, shearing, rotations, due to a com- 
plex balance between graphene's flexural modulus 
and interaction with Ir. 

The second scenario is sketched in Fig. [8b, which 
shows several 10 nm domains, commensurate, thus 
pinned to the Ir(lll) lattice, all having the same 
structure, but separated by domain walls, so-called 
discommensurations^i The strain field in these walls 
fixes their size: the 0.0093 A lattice parameter distri- 
bution FWHM suggests that one (or 1/3 or 2/3 if 
the domains are located on different terraces) Ir(lll) 
lattice parameter is accommodated in a few 100 
graphene unit cells, i.e. typically 100 nm. We also 
note that in this second scenario, the typical size of 
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FIG. 8. (Color online) Graphene deformations mapped 
onto a regular mesh with hexagons at its nodes, whose 
shapes represent amplified deformations (strain, rota- 
tion, shear), for two cases: ca. 60 nm domains, each 
having distinct lattice parameter (a), and a domain wall 
(ca. 100 nm) between two domains having the same lat- 
tice parameter (b). 



the graphene domains or domain walls is of the or- 
der of the average distance of ca. 100 nm between Ir 
step edges (Fig. [I]). This may not be fortuitous: Ir 
step edges, where the interaction between graphene 
and Ir is expected to be the strongest, could pin the 
graphene lattice, and drive the formation of com- 
mensurate graphene domains. Beyond a critical size, 
growing these commensurate domains would cost 
too high elastic energy with regard to the low inter- 
action with Ir, so that inhomogeneous deformation 
and shears/rotations resulting in incommensurabil- 
ity between steps are preferred. This would result 
in large discommensurations extending over the ter- 
races. The increase of inhomogeneous strain and 
rotations with decreasing temperature would result 
naturally from the increasing lattice parameter mis- 
match with respect to the growth temperature. This 
scenario would more easily explain than the first one 
the partial inheritance of the Ir TEC by graphene, 
as this could result from a complete inheritance in 
the small commensurate domains near step edges. 

Detecting the variations of lattice parameters dis- 
cussed above with STM is an extremely demand- 
ing task, and our attempts in this direction remain 
unsuccessful. In the fist scenario, this would re- 
quire detecting variations of ac of a few 0.001 A 
(the FWHM of the distribution of a c is 0.0093 A, 
see above), or of a few 0.01 A taking benefit of the 
ca. 10-fold amplification of strains by the moire. 
This is already below or at the verge of the reso- 
lution of the STM instrument which has been em- 
ployed whose resolution is limited to 0.004 A and 
0.04 A precisions for direct ac measurements and 
moire analysis respectively, 33 which are typical val- 
ues for high resolution STM. In addition, such varia- 
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tions should be detected between graphene domains 
distant by several 10 nm (first scenario), or contin- 
uously over ca. 100 nm-wide discommensurations 
(second scenario), which would require drift effects 
to be less than a few 0.001 A and 0.01 A for direct 
ac analysis and moire analysis respectively, another 
requirement which cannot be reached in most STM 
instruments. We note that micro-diffraction, such 
as possible in a low-energy electron microscope, has 
been used to detect ca. 0.002 A variations of ac 
in graphene/Ru(0001). 32 This technique does not 
however provide sufficient lateral resolution (set by 
the diameter of the aperture employed for defining 
the size of the micro-spot) for detecting the varia- 
tions in graphene/Ir(lll), which are characterized 
by length-scales in the range of several 10 nm. 



VII. CONCLUSION 

The carbon bond length in graphene on Ir(lll) 
is equal to that calculated for isolated graphene 
to within 0.2%. High resolution X-ray scattering 
measurements reveal that the structural parameters 
in very high quality graphene have some scatter, 
0.3% regarding lattice parameter and 0.5° regard- 
ing in-plane rotations. These small variations are 
either due to the presence of domains, with size ca. 
60 nm and different structure, most often incom- 



mensurate with Ir(lll), or to discommensurations, 
extending across ca. 100 nm, between graphene 
domains pinned to the substrate defects and com- 
mensurate with it. We anticipate that such struc- 
tural variations may also be encountered in other 
systems weakly interacting with their substrate, for 
instance graphene on Cu or dichalcogenates pre- 
pared by chemical vapor deposition. The struc- 
ture of graphene, average lattice parameter, rota- 
tions and strain, is found to depend sensitively on 
temperature, which could allow for epitaxial con- 
trol of its properties. Understanding the influence 
of the observed structural variations over particle or 
quasiparticle scattering, phonon modes, local charge 
transfers and chemical inhomogeneities, which steer 
many of graphene's properties, call for further inves- 
tigations. We anticipate that intercalation of species 
suppressing the non- vanishing graphene-metal inter- 
action, like alkali metals or hydrogen, could yield 
truly free-standing graphene where the structural in- 
homogeneities could be released. 
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